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Abstract 

We discuss the thermal (or gravitational) responses in topological superconductors and in topological phases in 
general. Such thermal responses (as well as electromagnetic responses for conserved charge) provide a definition of 
topological insulators and superconductors beyond the single-particle picture. In two-dimensional topological phases, 
the Streda formula for the electric Hall conductivity is generalized to the thermal Hall conductivity. Applying this 
formula to the Majorana surface states of three-dimensional topological superconductors predicts cross-correlated 
responses between the angular momentum and thermal polarization (entropy polarization). We also discuss a use 
of D-branes in string theory as a systematic tool to derive all such topological terms and topological responses. In 
particular, we relate the Z2 index of topological insulators introduced by Kane and Mele (and its generalization to 
other symmetry classes and to arbitrary dimensions) to the K-theory charge of non-BPS D-branes, and vice versa. We 
thus establish a link between the stability of non-BPS D-branes and the topological stability of topological insulators. 
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1. Introduction 



A vital role of topology in quantum transport phenomena in solids has been recognized as a driving force of 
a dissipationless current. The Berry curvature [1] in momentum and real spaces induces the velocity of Bloch 
electrons [2], resulting in various effects such as the anomalous Hall effect [3] and the spin Hall effect [4]. 
Even the heat current can be induced by the Berry curvature as observed in the anomalous Ncrnst effect [5] 
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and the thermal Hall effect [6]. These quantum transport phenomena of topological origin (or: topological 
currents) are regarded as a promising candidate for spintronics with low energy cost. 

In metallic systems the usual transport currents with dissipation are dominant over these topological 
currents, and hence it is rather difficult to identify the latter. A topological current, however, can be non- 
vanishing even in insulating systems - such systems on general ground are called topological insulator. This 
is so since it is the topological properties of electronic wavef unctions, rather than a band structure of energy 
levels in solids, that are responsible for generation of a topological current. 

The quantum Hall effect (QHE) is a canonical example of topological insulators characterized by a topo- 
logical (i.e., quantized) transport law [7,8,9]. Recently topological insulators realized by strong spin-orbit 
interactions in two and three dimensions have been discovered [10,11,12,13,14,15,16,17]. In particular, the 
three-dimensional (3d) topological insulator is characterized by the topological magnetoelectric effect [18,19] 
(the axion electrodynamics [20] ) . (See below for more details.) As a consequence of non-trivial electrical wave- 
functions in the bulk, these topological insulators support anomalous boundary (edge or surface) modes, 
whose gaplcss nature is topologically protected. 

Analogously, for superconductors (SCs) and supcrfluids, one can consider topological properties associated 
with the wavefunctions of fcrmionic quasiparticlcs; Within the BCS mean-field theory, the BCS (Bogoliubov- 
de Genne) Hamiltonian is a fermion bilinear in the Nambu spinor. The case where we have a quasi-particle 
gap everywhere in momentum space is an analogue of the band insulator. A topological SC is a SC with 
a full gap and topologically nontrivial quasiparticle wavefunctions. Canonical examples of topological SCs 
include, e.g., the 2d chiral p-wave SC [21,22]. 

In three dimensions, the B phase of superfluid 3 He was recently identified as a new topological SC (su- 
perfluid) [23,24,25,26,27]. A recent surface transverse acoustic impedance measurement reported in Refs. 
[28,29] revealed a signature of the surface Majorana fermion mode on the surface of 3 Hc-B. A copper-doped 
topological insulator (Cu x Bi2Se3) has been discussed as a candidate of a 3d topological SC [30,31,32]. 

The purpose of this article is to describe the response theory of topological phases with a special focus 
on their thermal, rather than electrical, transport. A motivation for the thermal response is that it is 
well-defined even for phases in which the electrical charge is not conserved, such as topological SCs, or 
topological phases in spin systems. While topological SCs (as well as topological insulators) can be defined in 
terms of a topological invariant built out of, within the BCS mean-field theory, quasiparticle wavefunctions, 
the response theory gives a physically measurable definition of topological phases, which can be largely 
insensitive to microscopic details - a lesson we have learned from the physics of the QHE, or the topological 
magnetoelectric effect in the 3d topological insulator. In particular, formulating the thermal response as a 
response to external gravitational field, we will derive a thermal analogue of the topological magnetoelectric 
effect, which uncovers an interesting cross-correlation between thermal and mechanical responses, in terms 
of the temperature gradient, and an applied angular velocity, respectively (Sec. 3). 

For the bulk of this article, we will focus on topological SCs in two and three dimensions, such as the 
2d chiral p-wave SC or 3 He B. However, we will also briefly describe how such response theory can be 
systematically constructed for a wider class of topological insulators and SCs in the "periodic table" [23,33,34] 
(Sec. 4). In fact, we will illustrate that responses in topological phases are closely related to quantum 
anomalies in field theories. Such characterization of topological phases in terms of anomalies is expected to 
incorporate arbitrary strong interactions as far as a bulk gap is not destroyed. 

Interestingly, such considerations lead to a natural link to topological objects in string theory - D-branes. 
We will demonstrate, by considering a particular configuration of D-branes, we can realize a field theory 
model of topological insulators and superconductors with desired discrete symmetries (Sec. 4). The stability 
criterion of D-branes against Tachyon condensation is in one-to-one correspondence with the classification of 
topological insulators and superconductors (i.e., the periodic table [23,33,34]). From condensed matter point 
of view, the D-brane construction can be thought of as a convenient tool that bridges K-thcory classification 
of fcrmionic Hamiltonians, and liner response theory. 
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2. electromagnetic response in the QHE and 3d topological insulator 



c5p(t,x) = 5f(t,x) = T ^-_ = —e^d t A v (3) 



2.1. the QHE in two dimensions 

Let us start this review by first illustrating the usefulness of the effective field theory of linear response, 
by taking the QHE as an example. The electromagnetic response of the quantum Hall fluid is described by 
the effective Chern-Simons action 

/off = S / dtd2xe ^ A ^ A ^ kez, (1) 

where is an external electromagnetic gauge field. The Chern-Simons action can be derived by considering 
a coupling of the (topological) insulator in question to the external (background) electromagnetic field A^, 
and then by integrating out fermions to derive the effective action for A^. 

The Chern-Simons action encodes all types of topological responses in the QHE. (i) The current Sj (t,x) 
induced by A^ is computed, for time-independent vector potential, as 

This is the QHE, J x = o-hE v , with the Hall conductivity an = ke 2 /(2irh). (ii) Similarly, one can compute 
the charge 6p(t, x) induced by A^ as: 

Slog _ e^fc 
5A (t,x) ~ 2nh 

where c is the speed of light. This means that a (solitonic) flux tube binds k units of charge, which is, in 
cylinder geometry, nothing but charge-pumping in Laughlin's thought experiment. For uniform magnetic 
field B z , the change in total electron charge is 

ID8 SQ e 5N e 

cbQ e = a H oB => cj^ = ce-^j = a H , (4) 

where N e is the total electron number. This is nothing but the Strcda formula [35,36]. (iii) The Chern-Simons 
term can be written, for static configurations of A^, as I e g = oh J dtcPx e ^ A^diAj = (er#/c) J dtd 2 x 4>B Z , 
where A$ =: <p/c. If we define (the change in) magnetization by 

6M* = ^ = (o- H /c)S<P ec S -^ = a H , (5) 

where /i = e(f> is the chemical potential \pN e = (n/e)eN e = 4>Q e \- 

2.2. topological insulator in three dimensions 

Similarly, for the 3d topological insulator, the electromagnetic response is encoded in the effective action 
[18,19,20] 

Qe 2 f 6e 2 f 

^ = 32^ J dtd " X = ^-j *** E B. (6) 

This "axion" term can be derived, similarly to the Chern-Simons action, by integrating out fermions in 
the presence of the background electromagnetic field. The Dirac quantization condition and time-reversal 
symmetry (TRS) restrict 9 to be quantized, 9 = 0,7r (mod 27r); 9 = for trivial insulators or vacuum 
whereas 9 = 7r inside topological insulators. 

As inferred from the axion term, the topological insulator features topological magnetoelectric (ME) effect 

SI e s 9 e 2 „ 9 a „ 

M = — — = E = E (7) 

SB n2hc it Air ' { ) 

SI eS _6 e 2 6 a 

SE n2hc 7r47r ' K ' 
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8ind = (4ji/c)y" H P 




Figure 1. The three-dimensional topological insulator in cylindrical geometry; (a) the induced magnetization follows the direction 
of the applied electric field; (b) the induced polarization follows the direction of the applied magnetic field. In both cases, 
time-reversal symmetry is broken on the surface, by putting, say, magnetic impurities, as indicated by red arrows. In (a), the 
magnetization is caused by the surface Hall current j H . In (b), the polarization is caused by an excess or a deficit charge on 
the top or bottom surface [see Eq. (4)]. 

or, quantized electromagnetic poloarizability, 

5Pi _ SMi _ 9_^_ 

where a = e 2 /(He) is the fine structure constant. Equations (8) and (7) are the 3d analogue of Eqs. (4) and 
(5), respectively. 

The magnetization M in Eq. (7), which follows the direction of the external electric field E, is generated 
by the surface QUE (Fig. la); when a topological insulator is in contact with a topologically trivial insulator 
(or simply vacuum), the 0-angle in the axion term jumps by tt at the interface. As the axion term is the 
total derivative of the Chern-Simons term, such interface (where TRS is weakly broken) is accompanied 
by the half quantized surface Hall current with an = ±e 2 /(2h) which generates a bulk magnetization M. 
Similarly, when an external magnetic field B is applied, according to Eq. (4), it induces an excess or a deficit 
of charge on the surfaces which are orthogonal to B - this is the source of a bulk electric polarization P in 
Eq. (8) (Fig. lb). 

3. gravitational response in topological superconductors 

3.1. gravitoelectromagnetism 

3.1.0.1. energy- gravity coupling We now develop the thermal response theory of topological phases, with 
an eye, in particular, on applications to topological SCs. Our strategy here is, based on Luttinger's idea 
[38], to follow as close as possible the linear response theory for the electric charge: When we study charge 
response, we consider a small external probe field for the charge density, Hj ~ J d d x 4>(x)p(x). The same 
strategy can be adopted since energy is conserved; as in the charge response, we can consider an external 
(fictitious) source term which couples to the energy (Hamiltonian) density s(x), Hj ~ J d d x 4> g (x)e(x)/v 2 , 
where we have introduced v which has the dimension of velocity to assign a proper dimension to (f> g (see 
below). 

The external source <pg can be thought of as a fictitious gravitational field in the presence of the Lorentz 
invariance, as it allows us to identify the energy as the mass, s(x) = m(x)v 2 , where m(x) is the mass density. 
For the gravitational theory of nature (not for Luttinger's fictitious gravity that is a mere device to develop 
a linear response theory) , v should be replaced by the speed of light c, whereas in condensed matter systems 
with emergent Lorentz symmetry, v is the "Fermi" velocity (see below for more discussion). The coupling 
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of the system to the potential (j> g can then be written as Hj ~ J d d x <f> g (x)m(x), where, from the analogy 
to the electromagnctism, mass can be thought of as a "charge" coupled to gravitational field. 

With Lorentz invariance, in more covariant language (in the Lagrangian language) the gravitational cou- 
pling with the (thermal) energy current is introduced as follows: The energy density e and the energy current 
j E are components of the energy- momentum tensor T^ v , e = T 00 and j E = vT° l . They are thus coupled to 
the variation of spacetime metric p M „ as —(1/2) J dtd d x \J—g T^ v 8g^ v in the Lagrangian, where g = detg^ 
- in a way analogous to the way the charge current j M couples with the external electromagnetic potential 
as - / dtd d xj^A^. 

3.1.0.2. gravitoelectric field A spatial gradient in energy density inevitably causes a temperature gradient, 
as one can infer from the thermodynamic equality dU = TdS as follows. (Here, U is the internal energy, 
S the entropy and T the temperature). Let us divide the total system into two subsystems (subsystems 1 
and 2). Equilibrium between the two is achieved when the total entropy is maximized dS = dS± + dS2 = 0. 
Since energy is conserved, dE 2 = —dEi, and hence dS\/dEi — dS 2 /dE 2 = 0. I.e, T\ = T 2 . Let us now turn 
on a gradient in the gravitational potential, so that the gravitational potential felt by Subsystem 1 and 2 
differs by 8cj) g . In this case, we have dE 2 = —dEi{\ + 8<j) g /v 2 ). This suggests the generation of a temperature 
difference T 2 = T\(l + 5<f> g /v 2 ). In other words, we can view the "electric" field E g associated with the 
gradient of <f> g , which we call "gravitoelectric field", as temperature gradient, 

E g := -V(j) g = v 2 T~ 1 VT. (10) 

3.1.0.3. gravitomagnetic field The analogy with electromagnetism can be further put forward - such for- 
malism is called gravitoelectromagnetism [39]. For our purposes to develop a linear response theory for 
thermal transport, we put an external gravity field which is infinitcsimally small. We can thus write the 
metric as g^ v = ?; M „ + where rj^ v is the metric of flat space-time. In the presence of matter, in the post 
Newtonian limit, we keep only hoo, ha, and hoi [which are of order 0(v~ 2 )], whereas hij = 0(t>~ 4 ) (i ^ j) 
are neglected [39]: 

ds 2 = v 2 (-l + 2<j) g /v 2 )dt 2 - -A g ■ dxdt + (1 + 2cp g /v 2 )dx ■ dx, (11) 

where we have introduced the gravitomagnetic potential, A g ^ := — (w 2 /2)/io;. For small h^ v , the Einstein 
equation can be linearized up to these non-zero components of h^ u kept in Eq. (11), and the resulting 
equation looks structurally identical to the Maxwell equation. (However, note that we are considering an 
external gravity field as a source which does not have its own dynamics) . As we have seen, the gravitoelectric 
field E g corresponds to temperature gradient. What does the gravitational analogue of magnetic field, the 
gravitomagnetic field, B g , correspond to? It turns out that B g can be understood as the angular velocity 
vector of rotating systems, 

B g := V x A g = vn. (12) 

For a system rotating with the frequency fl = £l z z, this can be understood by making a coordinate transfor- 
mation from the rest frame to the rotating frame in which the metric in the polar coordinates (t, r, ip, z) takes 
the form ds 2 ~ —v 2 dt 2 + 2fl z r 2 dLpdt + r 2 dip 2 + dr 2 + dz 2 . One can then read off, from the definition of the 
gravitoelectromagnetic field, the non-zero gravito gauge potential A? = (v/2)il z r. In Cartesian coordinates, 
A g = [v/2)Vt z z x x, and B g = V x A g = vfl z z. 

3.2. thermal Streda formula for 2d topological SCs 

We now use the formalism described above to study thermal transport of topological SCs in two and 
three dimensions. The goal here is to establish a gravitational analogue of the Streda formula (4-5) and the 
topological ME effect (7-8) [40]. 

Let us start with 2d topological fluid which does not have time-reversal symmetry. To draw a parallelism 
with the charge response, we introduce the moment of the thermal current a s My = Me ~ (/V e )M, 
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where M^ v = (x^T°" - x u T ^) , and M% = M 12 . (Similarly, the thermal current can be decomposed into 
3t = 3e ~ (M/ e )jj where j is the electric charge current.) Here the average has to be taken at finite 
temperature: (• • ■ ) = /( £ n)( n l ■ • • |n) where e n and \n) are eigenvalue and cigcnstate of the Hamiltonian 
H, and f(s n ) is the Fermi distribution function. 1 2 We can then obtain 



v 



DM, 



Ktt 



T 



H 2 dT 



(13) 



from the Kubo formula for the thermal Hall conductivity 3 [40]. 

To see the physical meaning of Me, note that the definition of is similar to the orbital angular 
momentum: U IV = (1/c) {x^T u0 - x v T^°) . Indeed, when there is a rclativistic invariance, the energy- 
momentum tensor can be symmetrized so that T^ v = T v ^, and thus Mg V = cL^ v . While this is not the 
case in condensed-matter systems in general, the (pseudo) Lorentz invariance can emerge in solids at low 
energies, such as in graphene or in topological superconductors where electrons or quasiparticles obey the 
Dirac or Majorana equation. In these systems the Fermi velocity v plays a role of c, and the tensor is 
related to the orbital angular momentum as Me = vL. Thus, k h = (v 2 /2)(dL z / dT). This is an analogue 
of the electromagnetic Streda formula (5). 

To derive "the other half" of the Streda formula [a thermal analogue of Eq. (4)], we note that the variation 
of free energy is given by dF = —SdT — L ■ dQ where S is the entropy. In terms of the gravitoelectric and 
gravitomagnetic fields, this can be written as dF = -(TS/v 2 )(v 2 T- 1 dT)-(L/v)-d(vn) = -(Q T /v 2 )(dc/) g )- 
(M e/v 2 ) ■ dB g , where Qt = TS is the thermal energy density, which couples to <fi g /v 2 . From the Maxwell 
relation, we thus obtain 



v 3 (dM%\ v 3 fdQq 



" H 2T\dcj> g ) Bz 2T\dB* g 



(14) 



This is the thermal analogue of the Streda formula for the charge Hall conductivity, in that Qt is the zeroth 
component of the energy current as eN e is in the charge current. With the pseudo-relativistic invariance, 
this can be written as 



dL z \ v 2 ( dS 



dT ) nss 2 \dn* 



(15) 



3.3. cross- correlated response of 3d topological SC 

The thermal Streda formula derived above for 2d topological fluid can be used to study the response of 3d 
topological SCs to the temperature gradient and the rotation; Given that the fcrmionic quasiparticles are 
fully gapped in the bulk, all (topological) transport phenomena, in the presence of a boundary (surface) , can 
essentially be discussed by looking at the surface transport. For simplicity, one can consider a cylindrical 
geometry (as in Fig. 1), apply the temperature gradient d z T or external rotation VL Z , and discuss the 
responses, which are mediated by the surface. (As in the electromagnetic responses, we weakly break TRS 
at the surface and hence the surface is gapped). Such thought experiments lead to, for the (induced) thermal 
polarization Pe defined by SQt = -V ■ Pg, and for the moment of the energy current Me, 



1. Below, we will consider the part of the thermal current carried by j E (and hence Me) specializing to the case of fi = 0. 
This is valid in discussing topological SCs. 

2. Our convention for <j>g, A g , and Me in this review differs from the one employed in Ref. [40]; there is a factor of two 
difference in defining 9 and A g [see Eq. (11)]. Accordingly, there is a similar factor of 1/2 difference in defining Me- These 
factors are due to the spin 2 nature of gravitons, and are convenient since the resulting linearized Einstein equation does not 
have such factors of two, and the similarity with the Maxwell equation is clearer. Also, in the definition of <p g , we have put an 
extra factor of v 2 . 

3. While in the bulk, in the absence of interactions, the Kubo formula gives us one of the most direct ways to get the Streda 
formula for a^j and kjj , the applicability of the Streda formulas is not limited to non-interacting systems, as they can be also 
derived assuming the presence of chiral edge states. See Ref. [40] for details. 



() 



M E 
Pe 



{Tk h /v 3 ) 
(Tk h /v 3 ) 




(16) 
(17) 
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The parallel between the electromagnetic and gravitational cases are obvious. Observe, however, that the 
gravitational response is not quantized as in the charge response, because of the presence of the velocity v. 
The situation is somewhat similar to the detection of the conformal anomaly (central charge) from specific 
heat in Id quantum systems at criticality; while the central charge C for a given quantum critical system 
in Id is a dimensionless universal parameter, it shows up in specific heat Cy together with the velocity v as 



3.3.0.4. possible experiments The thermal (gravitational) analogue of the Streda formula [Eq. (15)] and 
topological ME effect [Eq. (18)] can be tested experimentally. An external angular velocity fi 2 results in the 
change in temperature (in 2d topological SCs) and thermal polarization (in 3d topological SCs). If the heat 
capacity of the system is sufficiently small, these may not be difficult to measure. 

In three dimensions, the "dual" response, i.e., the response to the applied temperature gradient [Eq. (17)], 
can be detected by making use of the Einstcin-dc Hass effect. Let us assume a cylindrical 3d topological SC 
is suspended by a thin string and apply thermal gradient (as in Fig. la). This induces surface energy current 
with angular momentum L z , according to Eq. (17). By the conservation law of total angular momentum, it 
must be compensated by a mechanical angular momentum of the material, which can be directly measured 
in principle. 

Note that in both responses to il z and to temperature gradients, the part of the responses of our interest 
are contributions from the fermionic quasiparticles. They should be distinguished from the contributions 
from bosonic excitations such as vortices. If Vl z is larger than the critical angular velocity fi c i above which 
vortices are introduced in the bulk of the sample, an extra contribution to thermal polarization would be 
generated. 

4. anomaly ladder and D-branes 

4.1. integrating out fermions and chiral anomaly 

We have constructed the response theory of 3d topological SCs starting from the thermal Streda formula of 
2d topological SCs. We now discuss the response of topological SCs (and insulators) to external gravitational 
field from a field theoretical point of view [41,42]. Our goal here is to derive a gravitational analogue of the 
axion term in the electromagnetic response, and show that it is related to quantum anomaly (chiral anomaly). 
In this section, we will use natural units e — c = h = 1, and set the Fermi velocity v to be unity for simplicity. 

Let us work with an example; A canonical example of the 3d topological SC is the B phase of 3 He, which is 
described, in momentum space, by the following BdG Hamiltonian: H = (1/2) J c? 3 fc (fc)?-^(fc)\I / (fe) , where 
'l't(fc) = (cj fc , cj fc , c.j. _ k , _ fc ) is the Nambu spinor composed of fermionic creation/annihilation operators 
(c k /c s k ) of a 3 He atom with spin s and momentum k, and the kernel %{k) takes the following form: 



where M is the mass of a 3 He atom, ji is the chemical potential, and | A| is the amplitude of the pair potential. 
With the d-vector pointing parallel to momentum, d(k) — |A|fc, there is an isotropic gap everywhere on the 
3d fermi surface. The critical point at (x = separates topologically trivial (/i < 0) and non-trivial (/j, > 0) 
phases, which are characterized by an integral topological invariant of symmetry class Dili (the winding 
number) v = and v = 1, respectively [23]. We will henceforth set |A| = 1 and drop the 0(k 2 ) term in £(fc), 
£(fc) — > — jj, = m. With a suitable unitary transformation, the BdG Hamiltonian is written in terms of the 



C v = CttT/(3v). 




(19) 
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4x4 Dirac matrices {0(1=1,2,35 0} as % = —i Y^i=i a i&i + mP- Upon this linearizeation, fi can be thought of 
as a "Dirac mass term" . 

We introduce an external gravity field which couples to the Dirac Hamiltonian (Lagrangian) as, S[ip, tp, e] = 
J dtd 3 x y^gC, 

C = $e»vf (d M - ^u^cd) i> ~ m<iH>, (20) 

where fx, v, . . . = 0, 1, 2, 3 is the space-time index, and a,b, . . . = 0, 1, 2, 3 is the flat index; e a M is vielbein, 
and uifj, ab is a spin connection; £ a & = [y ai 72,] /(4i). The effective action 7 e g[m, e] for the gravitational field is 
then obtained from the fermionic path integral I e g[m, e] = — ihi J7)[ijj,ip] exp(zS'['0, i/j, e]). 

Below, we are interested in the topological term of the effective action, which is, in the imaginary time path 
integral, an imaginary part of the iS e g[m, e]. This part can be computed by making use of chiral anomaly 
as follows [43]. We first observe that the continuum Hamiltonian H enjoys a continuous chiral symmetry: 
we can flip the sign of mass, in a continuous fashion, by the following chiral rotation ip ip = e^ 75 / 2 -!//, 
j/j —> ip = -0'e^ 75 / 2 , under which the mass is rotated as 

m'(4>) = rae"^ 75 = m (cos 4> + 175 sin0) , (21) 

so that m'(4> = 0) = m and m'(4> = tt) = — m. Since m can continuously be rotated into — m, one would 
think, naively, 7 e ff[w] = 7 e fj[— m]. This naive expectation is, however, not true because of chiral anomaly. 
The chiral transformation which rotates m continuously costs the Jacobian J of the path integral measure, 

V[^,^}=JV[i>',^], (22) 

This chiral anomaly (the chiral Jacobian J) is responsible for the gravitational analogue of the axion term 
(the 0-term). The Jacobian J can be computed explicitly by the Fujikawa method [44] as 

1% = -\nJ= \\^^ J d t x y /=g<? de tir bcd R? aB f, (23) 

where i? = dw+wAw {R a b^u = Q^Zh — d v u>ab + [^^i^v] a b) is the Riemann curvature tensor; as in the 
electromagnetic response, the 0-angle is fixed to either 9 = or 9 = tt by time-reversal symmetry. The former 
corresponds to a topologically trivial state, and 9 = tt to a topologically non-trivial state. The theta term 
in the gravitational effective action (23) ("the gravitational instanton term") is an analogue of the axion 
term cx 9E ■ B in the electromagnetic effective action; There is an obvious structural parallelism between 
the electromagnetic and gravitational cases [45]. 

To make the connection with the existence of topologically protected surface modes, we note, when there 
are boundaries (say) in the a; 3 -direction at x 3 = L+ and at a; 3 = L_, the gravitational instanton term at 
the non-trivial time-reversal invariant value 9 = tt of the angle 0, can be written in terms of the gravitational 
Chern-Simons terms at the boundaries, = ^csU 3 =l + — ^csU 3 =l_7 where (i,j,k = 0,1,2) 

/cS = \ 'h li / d3xellktr (uifyuk + ^WiWjUk) , (24) 

with d = 1/2. This value of the coefficient of the gravitational Chern-Simons term is one-half of the canonical 
value (l/47r) x (c'/24) with d = 1/2. As discussed by Volovik [47] and by Read and Green [21] in the context 
of the 2d chiral p-wave SC, the coefficient of the gravitational Chern-Simons term is directly related to the 
thermal Hall conductivity, which in our case is carried by the topologically protected surface modes. 

4.2. anomaly ladder and periodic table 

All types of responses we have discussed so far (the QHE, the topological ME, and their gravitational 
(thermal) analogues) are related to quantum anomaly; the Chern-Simons terms in the QHE (both electro- 
magnetic and gravitational) is a manifestation of parity anomaly; the axion term can be derived, as we have 
seen above, from chiral anomaly. In fact, a wider class of topological insulators and SCs in the periodic table 
(Table 1) can be related to quantum anomalies. 
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Symmetry Class\d 





1 


2 


3 


4 


5 


6 


7 ■•• 


A 


z 





z 





Z 





z 


■•■ 


AIII 





z 





Z 





z 





z ■■■ 


AI 


z* 











2Z 





z 2 


Z 2 ■■■ 


BDI 


z 2 


z* 











21. 





z 2 •■■ 


D 


z 2 


z 2 


z^ 











2Z 


■■• 


Dili 





z 2 


z 2 


1.0 











2Z ■ ■ ■ 


All 


22 





z 2 


z 2 


Z* 








■ • • 


CII 





2Z 





z 2 


z 2 


z* 





■•• 


C 








2Z 





z 2 


z 2 


z<? 


■•■ 


CI 











2Z 





z 2 


z 2 


■L<> ... 



Table 1 

Periodic table of topological insulators and superconductors for the 10 symmetry classes in various spatial dimensions. Topo- 
logical insulators (superconductors) in the complex symmetry classes (A and AIII) are related to the chiral U(l) anomaly. 
The primary series of the topological insulators (superconductors) with an integer (Z) classification in the eight real symmetry 
classes are located on the diagonal in the table. In even space-time dimensions (odd space dimensions) they are predicted from 
the chiral anomaly in the presence of background gravity (Z^), and from the chiral anomaly in the presence of both background 
gravity and U(l) gauge field (Z*). The topological response of topological phases in odd space-time dimensions (TP and Z*) 
follows from their higher-dimensional ancestor (ifi and Z*), respectively. 



Let us consider topological insulators and SCs with an integer topological invariant located on the diagonal 
of the periodic table (Table 1), which we call the "primary series". (The Z2 topological phases which 
"descend" from the primary series can be called 1st and 2nd descendants [34].) For the primary series, the 
effective action for electromagnetic and gravitational response can be derived essentially by repeating the 
procedure discussed above; taking a Dirac representative in D = d+ 1 = even, coupling it to the background 
electromagnetic and gravitational fields, and then integrate out fermions; chiral anomaly in D = d+ 1 = even 
dimensions allows us to calculate the imaginary part of the action (i.e., the part which encodes topological 
part of the response). The result is summarized as follows: 

S\nZ = 2TTi ch(E)A(R) . (25) 



Mo 



D 



S In Z represents the change in the effective action under the infinitesimal chiral transformation; by inte- 
grating it one gets the topological action for the linear response. ch(E) denotes the Chern character of the 
vector bundle E, which is explicitly given by ch(£') = tr [e F ^ 27T ^] in terms of its field strength F. A(R) is 
the A-roof genus and takes the form 

i(i?) = l + I ^tr[i? 2 ] + ..- , (26) 

where R is the curvature two-form on the manifold AId- Since (25) measures the number of chiral fcrmion 
zero modes minus that of anti-chiral ones as follows from the chiral rotation, (25) is equivalent to the index 
theorem in mathematics. 

Topological terms for d = even can be derived from the topological terms in d + 1 dimensions considered 
above: They are all Chern-Simons type and obtained as a boundary contribution from a (d+ l)-dimensional 
topological term. 

One could check that the topological response derived from chiral anomaly is fully consistent with the 
periodic table; (i) the topological terms derived in this way preserve correct discrete symmetries for primary 
series; (ii) for symmetry classes which are realized as topological SC (i.e., no charge conservation), the 
anomaly polynomial predicts that there is no topological response for EM field; only gravitational (thermal) 
response exists. 

One could think of this as an alternative "derivation" of the periodic table (for the primary series) . 
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D(-l) DO Dl D2 D3 D4 D5 D6 D7 D8 D9 



type IIB Z OZOZOZOZOZ 
typo IIA ZOZOZOZOZO 
09" (typo I) Z 2 Z 2 Z Z Z 2 Z 2 Z 

Table 2 

Dp-brane charges from K-theory, classified by K(§ 9_p ), K — 1 (S 9— p ) and KO(§ 9_p ). A Z 2 charged Dp-brane with p even or p 
odd represents a non-BPS Dp-brane or a bound state of a Dp and an anti-Dp brane, respectively. 

4.3. topological phases and D-branes 

4.3.0.5. introduction We would now like to point out an interesting connection between topological phases 
in condensed matter and D-brancs. D-branes are topologically stable objects in string theory; At the level of 
classical (supcr)gravity, which is a low-energy effective field theory of string theory, D-brane can be visualized 
as a (p + l)-dimensional solitonic solution to the lOd Einstein equation (such D-branc is called Dp-brane). 
Besides such geometrical attribute in gravity theory, an important property of D-branes for our purposes 
is the fact that open string excitations on D-branes give rise to a gauge field theory. This dual nature of 
D-branes has been proved to be a useful tool to understand non-perturbative phenomena in gauge field 
theories, such as monopoles, Seiberg-Witten theory [48], etc. (Such dual nature of D-branes also plays a key 
role in AdS/CFT or holography. For our discussions in this article, however, we will not use AdS/CFT.) 

We will describe below, what D-branes are, how we can construct topological phases from D-branes, and 
why they give us an insight on the periodic table. In fact, the opposite is also true; topological phases 
deepen our understanding of the stability of D-branes. As we will argue below, the stability of the so-called 
"non-BPS" D-branes is directly related to the Z2 topological index of topological insulators. 

4.3.0.6. K-theory charge of D-branes As mentioned above, D-branes are a topologically stable object in 
sting theory; Dp-brane is a (p + l)-d solitonic object in lOd space-time of classical (super)gravity theory. 
What is behind the stability of D-branes is the fact that there is a "charge" associated to them. These charges 
are quantized, and hence they cannot change for a smooth deformation of field configurations; D-branes are 
thus stable. 

The charge of a stable D-brane can be either integer- or Z2-valued, depending on the types of D-branes 
and sting theory. This is summarized in Table 2. Observe that for type IIA and IIB string theory, "0" and 
"Z" appear in an alternating fashion. On the other hand, for type I sting theory, the way K-theory charges 
"0" , "Z" and "Z2" appear closely follows the Bott periodicity. Compare entries "type IIB" and "type IIA" in 
Table 2 with "complex" symmetry classes (A and AIII) in Table 1, and entries "type I" in Table 2 with eight 
"real" symmetry classes in Table 1. In fact, it was argued that D-brane charges are classified by K-theory 
[49,50,51]; one then cannot help speculating on a possible connection between D-branes and the periodic 
table of topological phases. 

For some cases where a D-brane has an integral charge, D-branes are an electric or a magnetic source 
of an Abelian p-form gauge field C( p ) = (l/p\)C^J i2 "' tip dx ltl dx lta ■ ■ ■ dXp p , the so-called "Ramond-Ramond" 
(RR) gauge field. An integral of the RR-gauge flux generated by a magnetic D-brane, JgdC^, measures 
the charge of the D-brane, where S is a hypersurface which encloses the D-brane. For D-branes with a Z2 
charge, while we do see K-theory charge exist, it is not possible to write it down as a quantized integral of 
the RR flux. 

4.3.0.7. D-brane configurations for topological phases For our purposes, we consider D-brane configurations 
which consist of two types of D-branes, a Dp-brane and a Dq-brane. They are located in parallel and do not 
intersect (Fig. 2) [52,53,54]. One can then ask what kind of (gauge) field theory is realized in such D-brane 
configuration. In string theory, quantum fields are realized as a vibration of a fundamental string. 

In our configuration, an open string can have its end points on D-branes. Let us first consider a string 
which has one end on a Dp-brane and the other on a Dg-brane. Analyzing the vibrations of such string, 
one finds a massive fermion in the string spectrum. The mass of the fermion is proportional to the distance 
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Dp-brane Dg-brane 




Figure 2. A typical D-brane configurations that realizes a topological insulator (superconductor). 

between the D-branes. For a relativistic fermion which is fully gapped by a mass term, following the periodic 
table of topological phases, one can discuss (under a suitable set of discrete symmetries) its topological 
stability against adiabatic deformation; one can "compute" its topological invariant or topological charge. 
On the other hand, for the D-brane configuration at hand, if it is stable, one can assign a topological (K- 
theory) charge to it, following the K-theory charge of D-branes. One can check these two types of topological 
charges agree precisely. This is the first implication of the correspondence between the periodic table and 
D-branes. 

The Dirac fermion is not the only quantum field realized in the D-brane configuration. Let us now turn 
our attention to a string which has its both ends on the Dp-brane. Such string vibration gives rise to a gauge 
field living on the D-brane. (Its gauge group depends on the type of sting theory and D-branes.) What is 
the dynamics of such gauge field on the D-brane? This is answered by the effective action of D-branes. The 
topological part (the so-called "Wess-Zumino term") of the effective action for a Dp-brane in a flat space 

A 

IS 

(27) 



Swz = V / Cfg+D A ch(E)A(TX) 



where integration (J x ) is over the (p + l)-d world-volume X of the Dp-brane. 5 In this action, C( g+1 ) is the 
background (q + l)-form RR gauge field (which is, in our situation, sourced by the Dg-brane). On the other 
hand, ch(E) and A(TX) depends on the field configuration on the Dp-brane; The U(l) gauge bundle on 
the brane is denoted by E; TX and NX are the tangent bundle and normal bundle. Plugging the RR-ficld 
generated by the Dq-branc, whose integral J dC^ q+ \) is quantized as it measures the topological charge of the 
Dg-brane, the Wess-Zumino term of the Dp-brane recovers precisely the response theory that we discussed 
in terms of anomaly. 

As an example, let us consider the case with p = q = 5; 6 this D-brane configuration realizes the 2d 
QHE. The vibration of an open sting stretching between the D-branes gives rise to a (2+l)d massive Dirac 
fermion, T~L = %p' (-'Ei=i y&idi + mcz)^ where tp is the two-component complex Dirac fermion field. On 
the other hand, the WZ action of Dp-brane is given by the Chern-Simons term: 



Swz 



1 



2(2tt) 2 
-1 

8^ 



X 



C (2 ) A F A F = 



1 



(2) 



AAF 



2(2tt) 2 
±m 



(7 ( a) A d(A A F) 



8dm 



YxZ 



A A dA. 



(28) 



4. For the discussion below, we will focus on the situation where the topological charge of D-branes is an integer, for 
presentational simplicity. 

5. This is obtained from a more general expression (A.l) by noting that in our situation TX ® NX is a trivial bundle, and 
hence A(TX (g NX) = A(TX) A A(NX) = 1. 

6. Since we can apply the T-duality equivalence which shifts the value of p and q by one [57], we can fix the values of p, say 
p = 5. 

7. This configuration has 6 ND directions. 
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Here, the integral is over the 6d world volume of the Dp-branc, which we split into the common directions 
of the Dp- and Dg-branes (Y) and the compliment thereof (Z), and we noted that the Dg-brane couples 
magnetically to the RR two form C(2), and hence the integral (2tt)~ 1 J dCp) = ±1/2 measures the RR- 
charge of the Dg-branc. (We have done a partial integration.) The sign ± in front of the Chern-Simons 
term corresponds to sgn Ax = ±1 (see Fig. 2). Similarly, from the WZ coupling <S\yz = J x C(2) A A(R) = 
T2 fx ^( 2 ) ^ ^ Iff ' we obtain the gravitational Chern-Simons term ~ w Aduj + |w 3 , as expected from the 
responses. 

Observe the structural parallelism between the topological terms in the response theory of topological 
phases, and the Wess-Zumino action of D-branes. In the former, the effective action looks, typically, as 
J e ff oc (topological invariant) x J d d+1 x (topological term in gauge theory), where "(topological invariant)" 
is the topological invariant of the topological phase in question, and "(topological term in gauge theory)" 
is the term of topological origin in gauge theories such as the Chern-Simons term, or the axion term. 
For example, J e ff oc Ch x J d 3 x ^ vp A^d v A p in the QHE, where Ch is the TKNN integer. In the D-brane 
construction of topological phases, the coefficient in front of the topological term is given by the integral of 
the RR-held, and measures the K-theory charge of the Dg-brane [55,56]. 

5. conclusion 

We have described, in the order of increasing spatial dimensions, from two, three, and to arbitrary dimen- 
sions, the theory of thermal response in topological phases. Emphasized is a close analogy to electromagnetic 
responses of topological insulators by adopting the language of gravitomagnetism, which revealed a cross- 
correlation between thermal and mechanical (rotational) responses. 

Let us close with discussion on effects of interactions. The fact that topological currents (either in elec- 
tromagnetic or thermal response) are related to anomalies in field theories suggests topological phases with 
topological currents should be stable against interactions. This can be seen simply by observing that a term 
of topological origin in the effective action, once exists, has its coefficient which is quantized (in some case in 
the presence of some discrete symmetry). Thus, small interactions should not destroy topological properties 
of a given topological phase. 8 

We have computed the theta term (both electromagnetic and gravitational ones) by making use of chiral 
anomaly (which arises as a Fujikawa Jacobian). Let us imagine repeating the same calculations in the presence 
of interactions such as ipipip (a coupling to a bosonic scalar field ip), or electromagnetic interactions. The 
Adler-Bardeen theorem says that the anomaly (the Fujikawa Jacobian) will not be altered even in the 
presence of such interactions. (Historically, this nonrenomalization of the chiral anomaly was important to 
predict the number of quark colors even at the time when the details of the strong interaction were not 
known. ) 

Inscnsitivity of anomalies to interactions can further be illustrated by the "anomaly matching condition" 
proposed by 't Hooft; an anomaly (i.e., topological current) can be computed in terms of cither the infra-red 
(IR) or ultra-violet (UV) degrees of freedom, and the results should be the same. For example, the degrees 
of freedom in solids at UV are of "free-fermion" type, while deep in the IR region, such description can be 
replaced by quasi-particles (such as excitons) that arise due to interactions. 

The discussion above naturally echoes in the D-brane construction of topological phases; the field theories 
realized by Dp-Dg systems come with, in addition to massive fermions we discussed, other fields and interac- 
tions among them. As we can understand more or less geometrically, the topological phases are nevertheless 
stable. 

While we have discussed the stability of non-interacting topological phases against interactions, it remains 
largely an open problem what is the nature of strongly interacting topological phases which arise solely 
because of interactions, if they exist at all beyond the fractional quantum Hall effect. The existence of a 



8. An alternative derivation of the electrical/thermal Streda formula mentioned in Footnote 3 also provides another link 
between topological currents and anomalies, suggesting the stability of topological currents against interactions. 
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topological current, however, should be a hallmark that we can use to characterize even for these putative 
"fractional topological insulators" . 
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Appendix A. a short course for D-branes 
A.l. what is D-brane? 

In string theory [57], we consider a string as a fundamental object instead of an elementary particle. The 
vibration of a closed string produces, among others, the gravitational field at low energies. On the other 
hand, the vibration of an open string, which looks topologically like an interval with two ends, produces, 
among others, gauge holds in the low-energy limit. 

A D-branc is defined as an object where open strings can end on it [57,59]; D-brane can be considered as 
an object giving a boundary condition to string vibrations - the name D-brane originates from "Dilchlct" 
brane. 

In terms of the lOd gravity theory, which is the low-energy limit of string theory, D-branes are a stable 
solitonic object which is quite massive. A Dp-brane extends in p-spatial directions and in time-direction; i.e, 
its world volume is (p+ l)-d. For example, a DO-brane and a Dl-brane look like a particle (called D-particle) 
and a string (called D-string), respectively. 

For our purposes, we will note the following properties (see Ref. [59]): (i) On a D-brane, an abelian gauge 
theory is realized. The fluctuations of the gauge field correspond to the fluctuations of the end points of open 
string. If we have N coincident Dp-branes, we obtain a non-abelian U(AT) gauge theory (ii) An intersection 
of two D-branes realizes massless fermions. (iii) D-branes are characterized by a K-theory charge. The last 
point will further be discussed below. 

A. 2. K-theory classification of stable D-branes 

Gravitons and gauge fields arc not the only fields that arise in supcrstring theory. In addition to them, 
the vibration of a closed string generates Abelian p-form gauge fields, the Ramond-Ramond (RR) gauge 
fields. (Here, we are focusing for a moment on a particular type of supcrstring theory, type IIA and IIB 
supcrstring theory (with or without orientifolds). Supcrsymmctric (or BPS) Dp-branes are charged under 
these gauge fields; they have the RR charges, and hence stable. A Dp-branc directly (electrically) couples 
to the (p + l)-form RR field Cr p+ iy An anti Dp-brane is defined to be the one with a negative RR charge. 
In type IIA superstring, p takes only even integer values i.e., p = 0, 2, 4, 6, 8 and in type IIB, p takes only 
odd ones p = —1, 1, 3, 5, 7, 9. 

Moreover, a Dp-brane couples to other RR g-forms with q < p in the presence of the gauge flux on the 
brane. This is clearly described by the following formula of the RR couplings of a Dp-brane [58]: 



where TX and NX are the tangent bundle (i.e., the vector bundle tangent to the D-brane world- volume 
X) and normal bundle (i.e., the vector bundle which is normal to TX). The U(l) gauge bundle on the 




(A.1) 
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branc is denoted by Er>. This formula (A.l) means that if there is a non-trivial gauge bundle on the Dp- 
brane, it is possible that there exist charges which correspond to lower dimensional D-brane charges. Such 
a configuration can be interpreted as a bound state of a Dp-brane and the lower dimensional D-branes. 

In fact, there are D-branes which do not have any RR charges and which arc nevertheless stable. They 
are not supcrsymmctric and are called non-BPS D-branes [60]. Also a system of a Dp- brane and an anti 
Dp-brane sometimes forms a stable bound state. Such a system is called a brane anti-brane system [61]. 
They typically exist in the presence of the special projection called the orientifold projection. These brane 
configurations exhaust all possible D-branes in string theory. The orientifold projection means that we 
require the invariance of string theory under the action Q q = I q ■ f2, defined by the product of the parity 
/ with respect to q spacial coordinates and the orientation reverse of the string world-sheet. The set of 
fixed points of I q is called the orientifold (9 — q) plane. 

Being a stable object in string theory, one could imagine what is protecting them to "decay" . It turns out 
that one can assign a K-theory charge to D-branes, which is the reason of stability. Indeed, the K-theory 
provides a very systematic classification of D-branes in string theory [50,51]. In mathematics, (topological) 
K-theory classifies vector bundles. More precisely, we start from a pair of two bundles (E\, E 2 ) on a manifold 
X and consider its difference. In other words, we introduce the identification 

{E 1 ®H,E 2 ®H)~(E 1 ,E 2 ). (A.2) 

This defines the K-group K(X). 

In string theory, this identification is naturally interpreted as follows. We start with a brane anti-brane 
system. The gauge bundles on the D-brane and the anti D-branc are regarded as E\ and E 2 . Typically a brane 
anti-brane system becomes unstable because the total RR charge is vanishing and it can pair-annihilate. In 
other words, there appears so-called a tachyon field in the open string between the brane and the anti-brane. 
The tachyon field has unstable potential energy and condenses like the Higgs effect, which makes the system 
decay into a lower-dimensional D-brane. This procedure is mathematically described by (A.2), which means 
that the charge is conserved under the tachyon condensation. If the gauge field configurations are the same 
i.e., Ei = E 2l then the brane and anti-brane are completely annihilated, and nothing remains after the 
tachyon condensation. If E\ 7^ E 2 , then eventually the system decays into a D-brane which corresponds to 
the difference between E\ and E 2 ■ We presented the K-theory classification of D-brancs in type IIA, type 
IIB and type I string theory, where we take X = 8 9 ~ p for a Dp-brane via a compactification procedure. 
Notice that type I string theory is defined to be the projection of type IIB string theory by an (SO type) 
orientifold 9-plane. For type IIA and type I, we need to use a different K-thcory called K^ 1 (X), which just 
shifts the dimension by one, and KO(X), which is the real valued version of K{X). 

In particular, if we ignore the torsion of K-group K(X), then it is known that K (X) is reduced to even- 
dimensional cohomology (Bi>oH 2l (X, Q). Indeed, this is explicitly given by the Chern character and this 
nicely matches with the RR coupling formula (A.l) [49]. On the other hand, the argument based on K- 
theory with torsion is quite general and includes the case where the D-branes do not have any RR charges 
as is so for the non-BPS D-brane. 
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